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Upregulation of TET1 and downregulation of 
APOBEC3A and APOBEC3C in the parietal cortex 
of psychotic patients 

E Dong, DP Gavin, Y Chen and J Davis 

Increasing evidence suggests that epigenetic dysfunction may account for the alteration of gene transcription present in 
neuropsychiatric disorders such as schizophrenia (SZ), bipolar disorder (BP) and autism. Here, we studied the expression of the 
ten-eleven translocation (TET) gene family and activation-induced deaminase/apolipoprotein B mRNA-editing enzymes (AID/ 
APOBEC) in the inferior parietal lobule (IPL) (BA39-40) and the cerebellum of psychotic (PSY) patients, depressed (DEP) patients 
and nonpsychiatric (CTR) subjects obtained from the Stanley Foundation Neuropathology Consortium Medical Research 
Institute. These two sets of enzymes have a critical role in the active DNA demethylation pathway. The results show that TET1, 
but not TET2 and TET3, mRNA and protein expression was increased (two- to threefold) in the IPL of the PSY patients compared 
with the CTR subjects. TET1 mRNA showed no change in the cerebellum. Consistent with the increase of TET1, the level of 
5-hydroxymethylcytosine (5hmC) was elevated in the IPL of PSY patients but not in the other groups. Moreover, higher 5hmC 
levels were detected at the glutamic acid decarboxylase67 (GAD67) promoter only in the PSY group. This increase was inversely 
related to the decrease of GAD67 mRNA expression. Of 11 DNA deaminases measured, APOBEC3A mRNA was significantly 
decreased in the PSY and DEP patients, while APOBEC3C was decreased only in PSY patients. The other APOBEC mRNA 
studied failed to change. Increased TET1 and decreased APOBEC3A and APOBEC3C found in this study highlight the possible 
role of altered DNA demethylation mechanisms in the pathophysiology of psychosis. 
Translational Psychiatry (2012) 2, e159; doi:1 0.1 038/tp.201 2.86; published online 11 September 2012 



Introduction 

Epigenetic modifications of chromatin, including methylation 
of genomic DNA, and acetylation or methylation of histones, 
have a critical role in orchestrating the expression of specific 
genes in postmitotic neurons. Among these epigenetic 
modifications, DNA (cytosine) methylation at specific gene 
promoter CpG-rich regions is an example of a stable 
epigenetic modification that ensures and maintains neuronal 
phenotype identity. 1 Increasing evidence suggests that 
epigenetic dysfunction of the genome, in particular alteration 
of histone and DNA methylation marks, may account for the 
misregulation of gene expression present in neuropsychiatric 
disorders, such as bipolar disorder (BP), schizophrenia (SZ) 
and autism. 2-9 Patients with psychosis, including BP and SZ 
patients, express an increase of DNA methyltransferase 
(DNMT) 10-13 and a wide variety of DNA methylation altera- 
tions. 6 We recently reported that the reduction of brain- 
derived neurotrophic factor (BDNF) expression in the parietal 
cortex of SZ patients was associated with an increase of 
5-methyl cytosine (5mC) at a BDNF promoter. 4 Furthermore, 
reelin, a trophic factor expressed selectively in telencephalic 
GABAergic neurons, is downregulated in SZ patients and this 
dowregulation is associated with increased methylation at 
specific CpG promoter sites. 2 ' 3 ' 14 However, inconsistent 
results with studies of reelin promoter methylation in relation 



to mental illness 6 ' 1 5-1 7 suggest that the DNA methylation 
status of specific promoter regions in postmitotic neurons is 
not stable, but in contrast is a rapidly reversible process 
maintained by the equilibrium between the DNA-methylating 
action of DNA methyltransferases and the demethylating 
action of an active DNA demethylation process 18 ' 19 (Figure 1). 
This theory has been supported by several independent and 
particularly interesting findings showing that 5mC mark on 
promoter CpG-rich regions of specific genes can be oxidized 
to form 5-hydroxymethylcytosine (5hmC) by ten-eleven 
translocation (TET) proteins in mammalian brain. 20-24 Recent 
evidence suggests that TET proteins, important contributors 
in DNA demethylation process, have essential roles in 
neurodevelopment and aging 25 and are involved in human 
cancers. 26 

Further, it has been proposed that 5hmC undergoes two 
successive processing steps: (a) a deamination step that is 
catalyzed by the AID (activation-induced deaminase)/APO- 
BEC (apoliporotein B mRNA-editing enzyme complex) family 
of cytosine deaminases, turning 5hmC into 5-hydroxylmethy- 
luracil (5hmU); and (b) the base excision repair pathway, in 
which 5hmU can be removed and replaced by unmethylated 
cytosine by a group of glycosylases such as MBD4 and 
TDG 4 27-29 Rapidly inducible growth arrest and DNA 
damage-inducible proteins 45 (GADD45) 30-32 are thought 
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Figure 1 Schematic representation of putative DNA methylation/demethylation 
pathways. AID/APOBEC, activation-induced deaminase/apolipoprotein B RNA- 
editing catalytic component; BER, base excision repair; C, cytosine; DNMT, DNA 
methyltransferase; 5hmC, 5-hydroxylmethylcytosine; 5hmU, 5-hydroxylmethylur- 
acil; 5mC, 5-methylcytosine; TDG, thymine DNA glycosylase; TET, ten-eleven 
translocation protein. 

to be master coordinators of this process by recruiting 
deaminases and glycosylases to 5mC and 5hmC-rich 
promoter regions 4 (Figure 1). 

Although numerous studies of DNA methylation have been 
carried out in the post-mortem brain of psychotic (PSY) 
patients, 2-4 ' 6 the different components of the DNA demethy- 
lation pathway, with the exceptions of GADD45 and MBD4 or 
TDG, 4 remain largely unknown. 

In this paper, we initiated an investigation on the expression 
of the TET gene family and AID/APOBEC deaminases in the 
inferior parietal lobule (IPL) and the cerebellum of a cohort of 
PSY patients, which includes BP patients with psychosis, and 
SZ patients, major depression (DEP) patients and nonpsy- 
chiatric (CTR) subjects obtained from the Stanley Foundation 
Neuropathology Consortium Medical Research Institute 
(Bethesda, MD, USA). 

Materials and methods 

Subjects. IPL (BA39-40) and cerebellar cortex were 
obtained from the Stanley Foundation Neuropathology 
Consortium. We selected IPL because this brain area is 
among the most seriously damaged cortical regions in SZ 
patients. 33 Imaging studies have reported decreased volume 
and differences in hemispheric asymmetry in patients 
compared with controls. 34 In addition, many of the IPL 
functions, such as sensory integration, body image, concept 
of self, visuospatial abilities and executive functions, have 
been found to be disturbed in SZ patients. 33 Also, GADD45 
proteins have been shown to be upregulated in this region of 
the cortex in PSY patients 4 

The cerebellum was included in this study because this 
brain area contains high levels of TET. 35 Cerebellar cortices 
were cut sagittally, lateral to the vermis and did not include the 
dentate nucleus. 

The demographics associated with each group of subjects 
are listed in Supplementary Table 1. The cohort of 'PSY 
patients' includes nine BP disorder patients with psychosis 
and 10 patients with SZ. The cohort of 'depressed (DEP) 
patients' includes 12 subjects with major depression. Normal 
control (CTR) subjects include 1 1 non-psychiatric subjects 
who had no history of psychiatric disorders. Among CTR and 
the patient cohorts there are no significant differences in terms 



of the mean age (F (2 , 41) = 0.657, P= 0.524) or the gender 
distribution (x 2 = 3.273, P=0.07). An one-way analysis of 
variance (ANOVA) did not reveal significant differences 
regarding brain pH (F (2 ,40) = 1 -572, P= 0.220), post-mortem 
interval (F (2 , 41) = 1.604, P=0.213) or RNA integrity 
(F (3j 38) = 2.526, P= 0.072) among the CTR subjects and the 
patient cohorts. Details regarding clinical characteristics and 
tissue handling for the post-mortem samples were described 
previously. 36 To confirm the results obtained in the IPL of 
Stanley Foundation Neuropathology Consortium in a different 
cortical area (BA9) and a different cohort (Harvard Brain 
Tissue Resource Center, Belmont, MA, USA), we measured 
TET1 mRNA in the samples of 7 CTRs and 1 2 PSY patients (6 
SZ and 6 BP). The samples were matched for age (CTR: 
51.0 ±10.5; BP: 47.3 ±16.6; SZ: 54.3 ±7.5 years), gender 
(CTR: 4M/3F; BP: 3M/3F; SZ: 3M/3F), post-mortem interval 
(CTR: 20.2 ±3.6; BP: 21 .5 ±3.8; SZ: 16.9±3.2h) and pH 
(CTR: 6.6 ±0.2; BP: 6.3 ±0.3; SZ: 6.8 ±0.4). 

Quantitative real-time polymerase chain reaction. Quan- 
titative real-time polymerase chain reaction was carried out 
using the Applied Biosystems Real-Time PCR System with a 
SYBR green master mix (Fermentas, Glen Burnie, MD, 
USA). Total RNA, isolated using TRIZOL reagent (Life 
Technologies, Grand Island, NY, USA), was further purified 
using the Qiagen RNeasy kit (Qiagen, Valencia, CA, USA). 
The RNA integrity number was measured with an Agilent 
2100 bioananlyzer (Agilent Technologies, Santa Clara, CA, 
USA). The primer sequences for the genes analyzed are 
summarized in Supplementary Table 2. Each sample was 
run in duplicate and repeated twice. For normalizing mRNA 
expression, several housekeeping genes (NSE, GAPDHanti 
fi-actin) were chosen as the internal control. For each 
housekeeping gene, we measured the gene stability ranking 
using the NofmFinder algorithm. 37 This procedure allows for 
the identification of the housekeeping gene best suited for 
normalization. Because each of the genes studied yielded 
similar results when normalized to either fi-actin, GAPDH or 
NSE, and because fi-actin has the highest housekeeping 
gene stability (NormFinder), we normalized our data to fi-actin. 

Detection and quantification of 5mC and 5hmC by 
immuno-dot-blot analysis. To measure 5mC and 5hmC 
expressed on genomic DNA, total DNA was isolated using a 
Puregene Core Kit A (Qiagen). The amount of DNA to be 
used in the dot-blot assay was determined by titration of the 
sample with various concentrations of standard DNA. In all, 
25 ng of DNA, which is within linear range, were chosen for 
each sample. DNA extracted from the samples and 12.5, 25, 
50 and 100pg of 5mC or 5hmC standard DNA (EF-01 1-001 
and EF-01 0-002; Diagenode, Denville, NJ, USA) were 
denatured by adding denaturing buffer (0.5 m NaOH and 
1.5 m NaCI) and incubating at 65 °C for 45min. After 
neutralization (0.5 m Tris-HCI, 1.5 m NaCI), the solution was 
loaded in duplicate into a well placed on a positively charged 
nylon membrane (Hybond; Amersham, Pittsburgh, PA, 
USA). The wells were then washed once with 2 x saturated 
sodium citrate. The membrane was removed and briefly 
washed with 2 x saline-sodium citrate buffer and dried. The 
crosslinking of the DNA with the membrane was carried out 
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with a UV Stratalinker (Stratagene, Santa Clara, CA, USA) 
on the Auto Cross Link setting. To measure DNA recovery, 
the membrane was stained with SYBR gold dissolved in 
Tris-borate-EDTA buffer for 10min. The membrane was 
scanned on a STORM 860 (Amersham). The fluorescence 
signal of DNA on each spot was measured and calculated 
based on the standard DNA of the same membrane. The 
membrane was washed with Tris-buffered saline buffer 
containing 5% milk for 1 h and incubated with a primary 
antibody of either 5mC or 5hmC (1:10000; Active Motif, 
Carlsbad, CA, USA) for 2h, and then incubated with the 
proper secondary antibody for 1 h. The immunoreactive dots 
were detected by the application of the ECL Plus chemilu- 
minescence kit (Amersham Pharmacia, Pittsburgh, PA, USA) 
on a STORM 860 fluorescent imager (Molecular Dynamics, 
Sunnyvale, CA, USA) and quantified with IMAGEQUANT 
software. The linear regression standard curves of 5mC and 
5hmC with ^ = 0.996 are shown in Supplementary Figure 1. 

Methylated DNA immunoprecipitation. Methylated DNA 
immunoprecipitation experiments were conducted using a 
previously published protocol. 4 5mC monoclonal mouse 
(Diagenode) and 5hmC (Active Motif) antibodies were used 
in methylated DNA immunoprecipitation. We decided to use 
the procedure because bisulfite and most enzyme-depen- 
dent methods are incapable of distinguishing 5mC from the 
approximately 14% of methylcytosines in the brain that are 
5hmC. 8 Recent studies demonstrate that these 5mC and 
5hmC may have very different functions and genomic 
locations. 38 ' 39 The percentage of methylated vs unmethy- 
lated promoter was calculated using the following equation: % 
(meDNA - IP/total input) = 2 A [(Ct(1 0% input) - 3.32) - 
Ct(meDNA-IP)]x100% (adopted from MagMeDIP kit 
instruction manual; Diagenode). 

Western blot analysis. The samples were homogenized 
directly in Laemmli buffer (100jil/10mg of tissue). The 
protein extract was separated by 4-12% Tris-glycine gel 
(Invitrogen, Grand Island, NY, USA) and transferred to 
nitrocellulose membranes (Invitrogen). The membrane was 
blocked for 1 h at room temperature with 5% non-fat milk in 
phosphate-buffered saline (10itim phosphate-buffered sal- 
ine, pH 7.4) and then reacted with 1:2000 anti-TET1 



polyclonal antibody (GeneTex, Irvine, CA, USA). The immuno- 
reactive bands were detected by using anti-rabbit IgG 
(Amersham) at 1:1000 for 1h followed by the application 
of the ECL Plus chemiluminescence western blotting kit 
(Amersham). The immunoreactive bands were detected by 
the application of the ECL Plus chemiluminescence kit 
(Amersham) on a STORM 860 fluorescent imager (Molecular 
Dynamics) and analyzed with IMAGEQUANT software. 

Statistical analysis. One-way ANOVA and analysis of 
covariance with Bonferroni post-hoc comparisons were 
performed using the mRNA content of the genes measured, 
5mC or 5hmC levels, TET1 protein content, RNA integrity 
number, age, gender, post-mortem interval, substance use 
and alcohol abuse, to evaluate differences among the three 
groups (patients with psychosis (PSY), depression (DEP) 
and normal controls (CTR)), using PASW v. 18 software 
(SPSS, Chicago, IL, USA). In all samples, P-values were 
two-tailed, and comparisons were considered to be statisti- 
cally significant for P<0.05 level with the Bonferroni multiple 
comparison correction. 

Results 

mRNA levels of TET1 , TET2 and TET3. First, we examined 
the mRNA levels of TETs (TET1 , TET2 and TET3) in the IPL 
of CTR subjects, and DEP and PSY patients. As shown in 
Table 1, ANOVA reveals that for TET1 mRNA, there is an 
overall significant difference among the three groups (F (2 , 
32) = 7.43, P= 0.002). We found an approximately twofold 
increase in the expression of TET1 mRNA in the PSY group 
compared with the CTR group (P= 0.001, Dunnett's correc- 
tion; P= 0.002 Bonferroni correction). The increase of TET1 
mRNA was not paralleled by a significant increase of TET2 
and TET3 mRNA. Furthermore, the DEP group failed to 
show statistically significant increase in TET1, TET2 and 
TET3 mRNA compared with the CTR group. The increase of 
TET1 mRNA expression in the PSY group seemed to be 
region-specific because in the cerebellum of the same 
patients (Table 1) we failed to find an increase of TET1 
mRNA compared with the CTR group (F (2 , 37) = 1.711, 
P=0.195). Table 1 also shows that the TET1 mRNA was 
expressed at a relatively higher level in the cerebellum than 



Table 1 Summary of the TET1 , TET2 and TET3 mRNA expression in the IPL (parietal cortex) and the cerebellum of CTR subjects and DEP and PSY patients 



Variable (mRNA) Patient cohort Overall ANOVA Multiple comparison 

1 CTR 2 DEP 3 PSY F d.f. P-value DEP vs PSY vs 



13 12 



Parietal cortex 






















TET1 


0.64 ±0.09 


1.02±0.18 


1.29±0.10 


7.656 
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0.002 






0.00.1 




TET2 


1.38±0.15 


2.0 ±0.24 


1.93 ±0.23 


1.786 


2 


NS 










TET3 


0.91 ±0.16 


1.23±0.14 


1.28 ±0.11 


2.187 


2 


NS 










Cerebellum 






















TET1 


3.45 ± 0.49 


2.95 ±0.44 


3.28 ±0.27 


1.843 


2 


NS 











Abbreviations: CTR, nonpsychiatric; DEP, depressed; d.f., degree of freedom; IPL, inferior parietal lobule; NS, nonsignificant; PSY, psychotic; TET, ten-eleven 
translocation protein. 

TET1 mRNA expression was significantly increased in the PSY group compared with the CTR group. Each value represents the mean ± s.e.m. of relative quantity 
(2" ACt x 10" 2 ) of TET mRNA vs p-actin mRNA 
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Figure 2 (a) Ten-eleven translocation (TET)1 protein expression in the 
nonpsychiatric (CTR), depressed (DEP) and psychotic (PSY) group. In the inset is 
the representative gel image of western blot of TET1 and p-actin. The TET1 protein 
expression was increased in the PSY group compared with the CTR group: *P< 0.01 
(one-way analysis of variance (ANOVA)). Each value represents mean ± s.e.m. of the 
optical density (OD) ratio of TET1 protein vs p-actin. (b) Levels of 5-methylcytosine 
(5mC) and 5-hydroxylmethylcytosine (5hmC) expression in total DNA isolated from 
the CTR, DEP and PSY groups. Each value represents the OD ratios of 5mC and 
5hmC immunostaining (calculated from their corresponding standard curves) and the 
DNA loaded on the membrane (calculated from standard curve). The 5hmC level was 
increased in the PSY group compared with the CTR group (*P<0.01), while 5mC 
levels failed to change. The data represent the mean ± s.e.m. 

in the parietal cortex. The observation that TET1 is higher in 
the cerebellum than IPL is consistent with findings from other 
research groups. 23 ' 40 The rational for combining SZ and BP 
disorder patients in a large PSY group was suggested by the 
observation that TET1 mRNA levels (SZ: 1.47 ±0.6; BP: 
1.31 ±0.41) and protein levels (SZ: 0.34 ±0.08; BP: 
0.30 ±0.07) were similar in SZ and BP disorder patients 
and higher than that in CTR subjects(mRNA: 0.64 ±0.09; 
protein: 0.07 ±0.01) (Figure 3c). 

To determine if elevated TET1 mRNA expression is 
associated with TET1 protein expression, we examined the 
TET1 protein level in the same samples. As shown in 
Figure 2a, the expression of TET1 protein was increased in 
the PSY and the DEP cohorts, but only the PSY group showed 
statistical significance compared with the CTR group (F (2 , 
39) = 5.89, P= 0.006). 

5mC and 5hmC levels in total DNA extracts. 5hmC is the 
oxidized product of 5mC catalyzed by the TET family of 
proteins. To examine if the levels of 5hmC in the DNA are 
affected by an increase in TET1 levels in the tissue, we 
measured 5mC and 5hmC levels by immunoblotting using 
specific 5mC and 5hmC antibodies (see Materials and 
methods for details). Figure 2b shows that 5mC levels in 
total DNA are substantially higher than the levels of 5hmC in 
each of the three groups; however, 5hmC levels are 
significantly higher in PSY patients when compared with 
CTR subjects, but not in the DEP group (F (2 , 36 ) = 5.47, 
P= 0.009; post hoc test: P=0.01). 5mC showed no 
significant differences among the three groups (F (2 , 
3 6 ) = 0.3, P= nonsignificant) (Figure 2b). 

Pearson's correlation analysis shows that the 5hmC levels 
were positively correlated with the level of TET1 mRNA in 
each of the three groups (Figure 3a) (P= 0.003). On the other 
hand, the levels of 5mC failed to show any correlation with 
the amount of TET1 mRNA (P=0.812) (Figure 3b). As 
most of the PSY patients had TET1 mRNA relative values 
over 0.8 and high levels of 5hmC (Figure 3c), these variables 
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Figure 3 Pearson's correlation analysis between relative ten-eleven translocation 
(TET)1 mRNA expression and 5-hydroxylmethylcytosine (5hmC) (a) or 5-methylcy- 
tosine (5mC) (b) levels of all groups, (c) Plots of TET1 mRNA and 5hmC levels from 
nonpsychiatric (CTR) subjects and psychotic (PSY) patients. The line represents 
median value. BP, bipolar disorder; DEP, depressed; SZ, schizophrenia. 

when considered together using binary regression analysis 
reinforce the ability of the two measurements to identify the 
PSY patients. The binary logistic regression takes into 
account both the TET1 mRNA and the 5hmC levels in 
classifying correctly all the CTR subjects and the PSY patients 
(x 2 = 25.5, d.f. = 2, P< 0.00001), except for one PSY patient 
misclassified as a control. Given the small number of samples, 
this finding must await replication, and it would take several 
large N samples to calculate statistically a more exact 
classification function. 

Higher 5hmC level at the GAD67 promoter in psychosis. 

A map of the genome-wide expression of 5hmC suggests 
that 5hmC is primarily expressed at CG dinucleotides located 
in the body (particularly exons) of the genes; however, 5hmC 
has also been found enriched in CpG islands located in the 
proximity of the promoter transcriptional star sites (TSS) 
where it may have an important role in transcription. 25 Here, 
we studied whether specific promoter regions of glutamic 
acid decarboxylase67 (GAD67), a gene whose transcription 
is downregulated in PSY subjects, 41-44 express an enrich- 
ment of 5mC and 5hmC at its CpG-rich TSS. We used 
methylated DNA immunoprecipitation to measure 5hmC and 
5mC at the GAD67 promoter using previously published 
procedures 4 As shown in Figure 4, we found significantly 
higher 5hmC levels at the GAD67 gene at -537 to -415 
and -145 to +21 in PSY vs CTR. In the same patients, 4 we 
reported an increase in 5hmC at the BDNF IXabcd gene at 
- 60 to +50 in PSY vs CTR, but not within BDNF exon IXd 
at +1185 to +1305. 5mC is also increased at the BDNF 
IXabcd promoter at - 60 to +50 in PSY vs CTR 4 However, 
we found no differences in 5mC at either the GAD67 
promoter site or within the BDNF IXd exon. As expected, 
there was a significant negative correlation between GAD67 
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Figure 4 Increased 5-hydroxymethylcytosine (5hmC) at the glutamic acid 
decarboxylase 67 (GAD67) promoter in psychotic (PSY) patients. Binding of 5- 
methylcytosine (5mC) (a) and 5hmC (b) to GAD67 promoter was measured by 
immunoprecipitation using specific 5hmC and 5mC antibodies with previously 
published procedures. 4 We find significantly increased 5hmC, but not 5mC, in PSY 
patients (□) compared with nonpsychiatric (CTR) subjects (□) at both regions of 
the GAD67 promoter examined. *P<0.05 vs CTR; **P<0.01 vs CTR. 

mRNA expression and GAD67 5hmC at -145 to +21 
relative to the putative TSS (Spearman = - 0.351 ; P= 0.04) 
and a trend for a negative correlation at -537 to -415 
(Spearman = -0.309; P=0.07). 

mRNA levels of AID/APOBEC deaminases. We next 
measured the gene expression of AID/APOBEC family 
genes because they are important cytidine deaminases that 
have an important role in the DNA demethylation path- 
way 24,45 The AID/APOBEC family is comprised of the AID, 
APOBEC1 , APOBEC2, APOBEC3 and APOBEC4 sub- 
groups 46 Quantitative real-time polymerase chain reaction 
failed to detect any reliable signal for AID, APOBEC1 and 
APOBEC2 in these samples, indicating a very low gene 
dosage in our IPL samples. In contrast, APOBEC3 and 
APOBEC4 mRNA are easily detectable. As shown in 
Table 2, ANOVA analysis revealed a significant decrease 
in the expression of APOBEC3A mRNA (F (2 , 2 5) = 15.38, 
P= 0.0001) in both the PSY and DEP groups (post hoc test: 
P= 0.0001 and 0.001, respectively). In addition, a dramatic 
mRNA downregulation was also observed in APOBEC3C in 
the PSY group compared with the CTR group (F (2 , 
23) = 5.71 8, P=0.01; post hoc test P=0.01), but not in the 
DEP group. Although APOBEC3B, APOBEC3F and APO- 
BEC3H levels in DEP and PSY groups, and APOBEC3D 
levels in DEP group tended to be lower than in CTR group, 
the decrease did not reach statistical significance. No 
changes in APOBEC4 mRNA expression were found in any 
of the three groups (Table 2). 

Potential confounding variables. To establish that the 
observed effects on the levels of TET1 mRNA, TET1 protein, 
APOBEC3A, APOBEC3C and 5hmC could not be attributed 
to potential confounding factors (age, brain pH, RNA integrity 
number, post-mortem interval), we analyzed the data using 
ANOVA and analysis of covariance. We found no significant 
differences among the three groups (the PSY or DEP patients 
and CTR subjects) with ANOVA. When the effect of these 
variables on the group differences of 5hmC, TET1 mRNA or 



TET1 protein was examined by entering them as covariates 
(analysis of covariance), we found the observed differences 
among the three groups were still significant. In addition, 
similar analyses of history of substance or alcohol abuse, 
antipsychotic, antidepressant or mood stabilizer treatment 
entered as covariant failed to influence group differences 
observed for 5hmC, TET1 mRNA and TET1 protein. 

Discussion 

The major finding presented in this paper is that TET1 (mRNA 
and protein) is markedly increased (two- to threefold) in the 
parietal cortex of PSY patients and this increase associates 
with an increase of 5hmC level at genomic DNA and with 
specific increases of 5hmC levels at GAD67 (Figure 4) and 
BDA/F 4 promoters in proximity of their TSS. Importantly, the 
increase of TET1 mRNA observed in the parietal cortex from 
the Stanley Foundation Neuropathology Consortium cohort 
was repeated in a small number of prefrontal cortex sections 
obtained from the HBTRC. We now are replicating these 
findings in a larger HBTRC cohort comparing 32 CTR and 42 
PSY patients. 

The increase may be specific to the cortex because the 
cerebellum of the same patients failed to show significant 
TET1 changes, even though the levels of TET1 in the 
cerebellum were 3-4 times higher than in the cortex. The 
reason why the cortex, but not the cerebellum, shows an 
increase in TET expression was not addressed in this study 
and remains only a matter of speculation. Other brain areas 
such as the hippocampus also require study. 

The increase of TET1 in the cortex of PSY patients is not 
associated with an increase of TET2 and TET3 mRNA and 
cannot be attributed apparently to confounding demographic 
variables, nor to the type and dose of antipsychotic used. 
Furthermore, the increase of TET1 and of 5hmC may not 
generalize to all major psychiatric disorders because it was 
absent in DEP patients. In this study, we analyzed the data by 
combining SZ and BP disorder patients in the PSY patient group, 
because the increase in TET1 mRNA was similar in the two 
groups (Figure 3) and because the pathogenesis and the 
symptomathology of these two diseases are related. This 
categorization is supported by recent findings showing 
that lower expression of GAD67 41-43 ' 47 and BDNF, 48 ' 49 
and overexpression of DNMT 11 and GADD45(3 4 are common 
biochemical alterations in both of these psychiatric disorders. 

Recent DNA demethylation studies suggest that TET1, 
by hydroxylating 5mC residues at CpG-rich gene bodies or 
at TSS of specific promoters, could have an important role 
in triggering neuronal activity-induced enzymatic removal of 
the repressive epigenetic DNA mark, 5mC, in mammalian 
brain. 22 

The role of TET1 may be to facilitate the removal of 5mC 
from CpG-rich gene regions via formation of the intermediate 
5hmC to favor the deamination catalyzed by AID/APOBEC 
family of cytidine deaminases and the removal of 5hmU by the 
base excision repair pathway. In an animal model, evidence is 
accumulating that TET1 and APOBEC are involved in 
neuronal activity-induced, region-specific, active DNA 
demethylation and subsequent activation of gene expression 
in adult mouse brain in vivo. 4,22,25,31 
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Table 2 APOBEC mRNA expressions 

Variable (mRNA) Relative APOBEC mRNA expression* Overall ANOVA Multiple comparison 

1 CTR 2DEP 3PSY F d.f. P-value DEP vs PSY vs 



13 12 



APOBEC3A 1.0 ±0.23 0.15±0.06 0.31 ±0.12 15.38 2 0.0001 0.001 — 0.0001 

APOBEC3B 1.0 ±0.05 0.74 ±0.07 0.42 ±0.02 1.26 2 NS — — — 

APOBEC3C 1.0 ±0.02 0.58 ±0.04 0.50 ±0.03 5.718 2 0.01 — — 0.01 

APOBEC3D 1.0 ±0.01 0.69 ±0.01 1.02 ±0.01 0.782 2 NS — — — 

APOBEC3F 1.0 ±0.12 0.80 ±0.08 0.88 ±0.07 0.173 2 NS — — — 

APOBEC3G 1.0 ±0.13 1.08 ±0.11 0.97 ±0.04 0.081 2 NS — — — 

APOBEC3H 1.0 ±0.05 0.86 ±0.04 0.68 ±0.03 0.561 2 NS — — — 

APOBEC4 1.0 ±0.10 1.02 ±0.17 1.09 ±0.11 0.198 2 NS — — — 



Abbreviations: AID, activation-induced deaminase; APOBEC, apolipoprotein B mRNA-editing enzyme; CTR, nonpsychiatric; DEP, depressed; d.f., degree of 
freedom; IPL, inferior parietal lobule; NS, nonsignificant; PSY, psychotic. 

Note that AID, APOBEC1 and APOBEC2 were not detected in the IPL of our cohort, but were considerably abundant in NT2 cell line (human neuronal progenitor cell), 
indicating that their expressions are extremely low in terminally differentiated neurons or glial cells. The relative basic expression levels (2" ACt x 10 " 3 vs p-actin) of 
APOBECs measured in the IPL of the CTR group are: 0.71 ±0.23 (APOBEC3A); 0.11 ±0.05 (APOBEC3B); 2.5 ±0.5 (APOBEC3C); 0.11 ±0.04 (APOBEC3D); 
0.42 ±0.97 (APOBEC3F); 0.55 ±0.16 (APOBEC3G); 0.15 ±0.07 (APOBEC3H) and 1.6 ±0.36 (APOBEC4), indicating that the APOBEC3 subgroup genes are 
expressed heterogeneously in the IPL. 

a The expression of APOBEC3A was decreased both in the PSY and the DEP groups compared with the CTR group, while the expression of APOBEC3C mRNA was 
decreased only in the PSY group. Each value represents the mean ± s.e.m. of relative quantity (2" AACt ) of APOBECs vs p-actin. 



Although the TET1 -dependent active DNA demethylation 
and gene expression upregulation may occur in the brain 
under physiological conditions, 22 in the IPL of our PSY patient 
cohort we found an increase of TET1 that positively correlates 
with an increase of 5hmC at genomic DNA (Figure 3a) and at 
gene promoters such as GAD67 (Figure 4) and BDNF, which 
have been consistently associated with downregulation in the 
expression of these target genes 41 ~ 43 ' 47-49 A possible 
explanation for the unexpected finding is that an increase of 
TET1 associated with gene expression downregulation is 
associated with a downregulation of some of the main 
APOBEC-deaminating enzymes. Our studies on APOBEC 
show that the most abundant isoforms (3A and 3C) of these 
cytidine-deaminating enzymes may be reduced significantly 
in psychosis preventing the active demethylation of the gene 
promoters, resulting in an accumulation of 5hmC, which by 
recruiting MBD3-NURD complex appears to be a DNA mark 
that suppresses transcription. 50,51 

Up to now there have been 13 human cytidine deaminases 
identified 46 In our parietal cortex cohorts, AID, APOBEC1 and 
APOBEC2 were virtually undetectable, and APOBEC3B, 
APOBEC3D, APOBEC3F, APOBEC3G and APOBEC3H were 
expressed at a very low level. Only APOBEC3A, APOBEC3C 
and APOBEC4 were detectable at a significant level. The 
expression of APOBEC3A and APOBEC3C were significantly 
decreased in the cortex of PSY patients (Table 2). 

In a previous study, we found no changes in the expression 
of the base excision repair pathway enzymes in PSY patients 
(that is, TDG, XPG, MBD4). 4 Hence, the decrease of 
APOBEC cytidine deaminase enzymes in the cortex of PSY 
patients may be the limiting factor that reduces the rate of 
conversion of 5hmC into 5hmU preventing demethylation of 
5hmC genes and inducing the gene transcriptional repression 
that may underlie the pathophysiology of major psychosis. 

An alternative explanation for the role of TET1 in transcrip- 
tional repression is that the repressive function of TET1 is 
independent of its catalytic activity. In fact, it has been found 
that TET1, which contains CXXC DNA binding sites, 52 ' 53 
recruits the polycomb repressive complex 2 and SIN3A co- 



mpressor proteins at target genes, suggesting that this 
repressive protein and SIN3A may have an important function 
in TET1 -mediated gene repression. 54 Our analyses of history 
of substance or alcohol abuse, antipsychotic, antidepressant 
or mood stabilizer treatment, entered as covariants in analysis 
of covariance, failed to reveal an influence on group differ- 
ences in 5hmC, TET1 mRNA and TET1 protein, but our con- 
clusion of no effect of these variables requires further validation. 
Although our statistical analysis with medications as a covariant 
failed to show an influence of drug treatment on TET and 5hmC 
expression, medication factors on gene expression in post- 
mortem studies cannot be excluded. In fact, the medication 
history in the demographic records lacks detail and precision. 
In addition, there is evidence that different antipsychotics 
have different effects on epigenetic mechanisms 41 

In summary, the increase of TET1 and the decreases in 
APOBAC3A and APOBAC3C found in this study together with 
previous studies showing an increase in DNMT and GADD45|3 
in the cortex of PSY patients support the hypothesis that an 
epigenetic misregulation of gene expression may be operative 
in the pathogenesis of PSY disorders. Clearly, to test the 
hypothesis there is a need for a more detailed investigation of 
the binding of TET1 , APOBEC3A and APOBAC3C, and the 
presence of 5hmC across the genome. In addition to 
consolidating the findings from human samples and to study 
drugs that may normalize the altered epigenetic mechanisms, 
we need to develop a rodent model that will allow us to 
investigate changes in TET1 , APOBEC, DNMT and GADD45(3 
similar to the changes detected in the brain of PSY patients. 
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